Detailed Stability Analysis of Electroweak Strings 
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ABSTRACT: 

We give a detailed stability analysis of the Z-string in the standard electroweak 
model. We identify the mode that determines the stability of the string and 
numerically map the region of parameter space where the string is stable. For 
sin^ Oyy = 0.23, we find that the strings are unstable for a Higgs mass larger than 
23GeV. Given the latest constraints on the Higgs mass from LEP, this shows that, if 
the standard electroweak model is realized in Nature, the existing vortex solutions 
are unstable. 



Introduction 



Recent studies have shown that the stabihty of topological defects may persist 
when they are embedded in theories where they are not topologically stable. Such 
"embedded defects" are exact solutions of the equations of motion in the theories 
where they are embedded and can be dynamically stable. 

A typical example of such a defect is the semilocal string [1,2], an embedding 
of the Nielsen-Olesen vortex in a "semilocal" model with SU {2) gi^iig^i x t/(l)/oea/ ^ 
U {X) global symmetry breaking. The semilocal string has been shown to be stable for 
a finite parameter sector [3,4] even though the vacuum manifold in the semilocal 
model is S'^ . In fact, since 7ri(5''^) = 1, the stability of the semilocal string is 
dynamical rather than topological. The crucial reason for this stability is that for 
a certain parameter region, the increase in gradient energy necessary for the string 
to decay is more than the corresponding decrease in potential energy. Therefore, 
for that parameter region, the decay would increase the total energy and is not 
favored energetically. 

The obvious generalization of the semilocal model is the electroweak model 
in which the SU{2) symmetry becomes gauged: SU{2)i x U{1)y U{l)em- It 
is possible to show that the Nielsen-Olesen vortex may be embedded in the elec- 
troweak model i.e. the electroweak equations of motion with a vortex-like ansatz 
reduce to the Nielsen-Olesen equations [5,6]. Two such embedded vortex solutions 
- the so-called r— and Z— strings - are known [7]. Here we shall only consider the 
Z— string as the r— string is expected to be unstable for all values of the parameters 
[7]. 

The phenomenological successes of the standard electroweak model [8] make 
the question of stability of the electroweak string a very interesting and important 



one. A stable electroweak string would for the first time in particle physics pro- 
vide a macroscopic, stable coherent state at low enough energies to be accessible 
to particle physics experiments. It is this question of stability [9] that we are ad- 
dressing in detail in this paper. In particular, we construct a map of the parameter 
space of the standard electroweak model showing the range of parameters where 
the Z— string is stable and where it is unstable. Here we give the details of the 
calculation as well as the physical reasoning for the simplifications that occur in 
this originally highly complicated problem. A brief exposition of the main results 
we derive here may be found in Ref.[15]. 

The structure of the paper is the following: in the next section we give a 
review of the electroweak string showing that the electroweak equations of motion 
reduce to the Nielsen- Olesen equations for a particular vortex ansatz. In section 
3 we show how can the stability problem of the electroweak string be reduced to 
the eigenvalue problem of a single Schroedinger-like eigenvalue problem. This is a 
non-trivial simplification since the initial system of coupled perturbations involves 
twenty coupled degrees of freedom which after tedious manipulations not only 
decouple but also reduce to a single eigenvalue equation. Finally in section 4 we 
solve this eigenvalue problem and construct a map showing the parameter sector 
corresponding to stability. 
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Review of Electroweak Strings 

We consider static bosonic field configurations in the Weinberg-Salam model: there 
is no time dependence and we choose a gauge where the zero components of the 
gauge fields are set to zero. The energy functional is given by 



/" 1 X 



(2.1) 



The Weinberg angle is given by tan 6w = g' /g. The masses of the W-boson, 
the Z-boson and Higgs boson are, respectively, 

Mw^\gV, Mz^^aT^, Mh ^V2X7^, (2.2) 

where 



« ^ + 9'' (2.3) 
The time-independent field equations are 

DjF^^j = - li^(0t r^Di<p - {Di<j>)^ T» (2.4) 

djfij = -W (fDict> - {Di(t>)U) (2.5) 

DiDi4> = 2X{4>U - \ri^)4> , (2.6) 



The symbols are in the standard notation defined in Ref. 10. In addition, we recall 
the usual mixing formula: 

Z^' = cos%W^^^ - sinOwB^" , A^" = sinOwW^^ + cosOwB^" , (2.7) 



The vortex solution extremising (2.1) is given by [5,6]: 

]nO = ' ' eg 



W^"^ = = 1^^ = A^', Z^' = [A^']J ^ivo(0. 



(2.8) 



where, the coordinates r and 9 are polar coordinates in the xy— plane. The integer 
m is the winding number of the vortex and, here, we shall restrict ourselves to the 
case m = 1. The subscript NO on the functions / and A^^ means that they are 
identical to the corresponding functions found by Nielsen and Olesen [11] for the 
usual Abelian-Higgs string. On substituting eq. (2.8) into the equations of motion 
they reduce to 

- 7 + « (i - f ") /' = (2-10) 

where primes denote differentiation with respect to r and the subscripts NO have 
been dropped for convenience. These are solved together with the boundary con- 
ditions: 

/(O) = = ^(0), /(oo) = ^, ^(oo) = - (2.11) 

V2 Oi 

The string solutions resulting from these equations have been studied previously 
by several authors in a lot of detail. A sample of these papers may be found in the 
collection of Ref. 12. 

At this point, it is useful to note the symmetries of the string configuration. 
Firstly it is axially symmetric i.e. it is invariant under the action of the symmetry 
operator generated by the generalised angular momentum operator 

K, = L, + Sz + h . (2.12) 

Lz and Sz are the usual orbital and spin pieces respectively of the spatial angular 
momentum operator. Explicitly 

d 

Lz = 1 {Szo)j = -iesjkSk 1, (2.13) 
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where a is any vector field and 1 is the 2x2 unit matrix. . Note that Sz annihi- 
lates the scalar Higgs field. Iz is composed of a U{1) generator, Y and an SU{2) 
generator, 



Y and T"^ act on the Higgs field on the left. Y annihilates the gauge field and T'' 
acts via a commutator bracket. 

The configuration has two further symmetries. It is invariant under the com- 
bination of reflection in the x-axis and complex conjugation. Also it is invariant 
under the action of the global U (2) gauge tranformation given by 



Our eventual expansion of the perturbations will be in Fourier modes but we shall 
point out some connections between these and the eigenfunctions of the above 
commuting symmetry operators. 



liY-T^). 



(2.14) 




(2.15) 
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Stability of Electroweak Strings 

The vortex solution given in eq. (2.8) is not topologically stable. This means 
that any field configuration can be continuously deformed to the vacuum. Hence 
we are investigating the metastability of the vortex solution i.e. whether it is 
a local maximum or minimum in configuration space. We consider infinitesimal 
perturbations of the vortex configuration and ask if the variation in the energy is 
positive or negative. 

Let us write 

Z^' = Zl^^ + SZl' (3.2) 

= diag{-cos2ew, 1), (3.3) 

and, 

dj = (djl + i-aT^Zj) . (3.4) 

The perturbations can depend on the 2;— coordinate and the 2;— components of 
the vector fields can be non-zero also. However, since the vortex solution has trans- 
lational invariance along the 2;— direction, it is easy to see that the 2;— dependence 
in the perturbations can be ignored and the 2;— components of the gauge fields can 
be set to zero. This follows from (2.1) where the relevant 2— dependent terms in 
the integrand are: 

iGhGh + ^FBisFBis + {Dscj>)HD^cl>) (3.5) 

This contribution to the energy is strictly non- negative and is minimized (that is, 
made to vanish) by setting the 2;— components of the gauge fields to zero and also 
considering the perturbations to be independent of the 2— coordinate. For this 
reason, we shall drop all reference to the 2;— coordinate in the calculations below 



and it will be understood that the energy is actually the energy per unit length of 
the string. 

Now we calculate the energy of the perturbed configuration discarding terms 
of cubic and higher order in the infinitesimal perturbations. We find, 



E = {Eno + SEno) + Ei + Ec + Ew 



(3.6) 



where, E]\jq is the energy of the Nielsen-Olesen string and SEjsjq is the energy 
variation due to the perturbations 02 cind SZ^^. The variation Ei is due to the 
perturbation 0i in the upper component of the Higgs field: 



El 



|(ij0l|2 + 2A(/2- 772/2)1011 



where, 



dj = dj — i—cos{29'\Y)Zj . 
The contribution from the (p and W"^ [13] interaction is: 



(3.7) 



(3.8) 



Ec = cosOw I d^xJfW^ 



1 



Jf = -ia 
^ 2 



and the energy in the W"" and A bosons is 



(3.9) 
(3.10) 



E-^Y = I d X 



X #2 . y X Z + -|V X #1 + X 



(3.11) 



+ -|V X #2 + 7Z X + \g^f{W°'f + -(V X If 



where, 7 = gcosOy^. It may be remarked that the / and Z fields in eqs. (3.7)-(3.11) 
are the unperturbed fields of the string since we are only keeping up to quadratic 
terms in the infinitesimal quantities. 

Firstly we note that the perturbations of the fields that make up the string 
do not couple to the other available perturbations, i.e. the perturbations in the 
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fields / and v only occur inside the variation SEjqQ. We can understand this as 
follows: the perturbation of the string solution has t/ = 1, where U is given by eq. 
(2.15)., whereas the other perturbations have U = —1. Now, since we know that 
the Nielsen- Olesen string with unit winding number is stable to perturbations for 
any values of the parameters then necessarily, SEj^fQ > and the perturbations 
02 and SZf^ cannot destabilize the vortex. Then, we are justified in ignoring these 

— * 

perturbations and setting 5Ej^q = 0. Also we note that the A boson only appears 
in the last term of eq. (3.11) and obviously makes a positive contribution so we 

— * 

can set A to zero. 

We now consider an expansion of the remaining perturbations in Fourier modes. 
This gives. 



01 = X{r)e 



imO 



(3.12) 



for the m^^ mode where m is any integer. For the gauge fields we have. 



{fl{f)cos{n9) + /ism(n^)} H — [—hisin{n9) + hicos(n9)^ bq 



(3.13) 



{— /2(r)sm(n6') + f2Cos{n6) } H — {^h2Cos{n6) + h2sin{n6) } cq 

for the n^^ mode where n is a non-negative integer. Inserting the expressions for 
the wf^ mode of 0i and the n^^ mode of the W°' fields in the energy functional 
gives: 

.2 



(3.14) 



Ei = 2t: j drr x'V |^ (m + |cos2%^)^ + 2A (^f - | 



X 



(3.15) 



Ec = 5±n,l-mT^occosew / drr -(/%'- x/)(/2 T /l) 

jipn + 2 - |(1 - cos2ew)v^ (/i2 ± hi) 
for n 7^ 0. In the case when n = 0, the expression for Ec is equal to the above 

9 



(3.16) 



expression multiplied by a factor of 2 and with /i and h2 set equal to zero. Finally, 

ypif2hi - fih2)v' 

+ ^ \ ~nfi + hi' -^vf2f + ^ \nf2 + 112 +^vfif (3- 

+ Ig^fVifl^ + h^) + + h2^]] + (/a ^ /a, h-a ^ /ia) • 
for n 7^ 0. In the case when n = 0, Et^ is given by (3.17) multiplied by a factor of 
2 and with fi and h2 set equal to zero. 

First let us inspect Ei . Here the negative contributions can come from the term 
proportional to . The coefficient of is composed of two terms: the second 
term comes from comes from the potential part and is negative while the first term 
comes from the kinetic part and is always positive and is smallest when m = 0, at 
least in the region near the center of the string where an instability is most likely 
to develop. That is, the m = mode is the "most dangerous" mode. 

Next we inspect E^/-. Here the analysis is less obvious. Yet one can see that 
the only term that can be negative is the first term that arises from the term 
X Z in eq. (3.11). For this to contribute at the center of the string, 
the vector fields must not vanish there. But the only mode that need not vanish at 
the center is the n = 1 mode; all other modes have to vanish at the center if they 
are to be single-valued and finite. Therefore the only mode that can give negative 
contributions at the center of the string is the n = 1 mode. Hence, we restrict 
ourselves to considering the m = 0, n = 1 mode. 

There is also a generally accepted idea which leads us to expect the least stable 
mode to be the m = 0, n = 1 mode. As noted, the original configuration has 
Kz = and U = 1 . The m = mode is also Kz = with U = —1. Turning our 
attention to the gauge fields since 

Kz e^^-^Ti = (n ± 1) e^^'^Ti , (3. 
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where t± = ri ± zt2, we see that taking the combinations /i — /2 and hi + h2 
we get Kz = and U = —1. Intuitively we expect the least stable mode to have 
Kz = since as \Kz\ increases one gets an increasing centrifugal barrier. The other 
combinations /i + /2 and hi — h2 are a superposition of Kz = ±2 and we will see 
they decouple and drop out of our analysis. 

Also note that the barred and unbarred perturbations decouple. ( Under the 
combined operation of reflection in the x-axis and complex conjugation the barred 
and unbarred variables have opposite sign.) Furthermore, the stability problem in 
the barred variables is contained within the problem of the unbarred variables (all 
we need to do is to set 0i = 0). Therefore, it is sufficient to consider only the 
unbarred functions. 

We will now systematically simplify the expression for the energy variation. 
After a lot of algebra, we obtain the first step: 



6E 
2^ 



dr r 



dr 



+ acosOw J dr r ^ — {fx ~ / X) - ^(1 - oisin 



+ T+(F_,x,^+) + T_(F+,^_) 
where. 



(3.19) 



r+(F_,x,e+) = 



dr 
2^ 



P+F- + 



{ii-iv)e++iv'^+} ^2(/x'-/'x)^2 



+ 



acosO^^r^ 



+ 



(3.20) 
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P± = (1^7^^) + 



2 , g'r'f 



= ^cosHOwv'^ + 2A ( - 



A±{r) = 
V 



/2±/l 
2 

h2 ± /ii 
2 



2P±(r) 

2 'A_M£ 
2 / 2P+ 



1 r gV^ff 
r dr \ 2P+ 



(3.21) 

(3.22) 
(3.23) 
(3.24) 

(3.25) 
(3.26) 



2i-2 



S±ir) 



± r— 
dr 
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v' (1 =F 7v) 
r P±{r) 



2 P±(r) 

As expected the problem in x, F_ and ^+ has decoupled from the problem in 
F+ and Furthermore the only terms containing F- and are T_ and 
respectively and since these are whole squares they can be set to zero. This then 
leaves us with a problem in x and ^-|- and a problem in just 

We first discuss the ^_ problem. We conjecture that the relevant potential, 5'- 
is positive for any values of the parameters (3 and cos6'p;j^ (where (3 — SX/g"^). We 
motivate this by considering the asymptotic behaviour as r — > and as r — > oo: 
in both these limits S- is always positive. Our conjecture is therefore reasonable 
since the function is basically exponential, modulated by polynomials. (We have 
checked it numerically for many pairs of parameters.) 

Hence our analysis is reduced to considering perturbations in ^-|- and x alone. 
We express the change in the energy in the form of an eigenvalue problem: 



(3.27) 



6E[x,^+] = 27T J drr(x,e+)0 



X 



(3.28) 



where, O is a 2 x 2 matrix differential operator. 
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It is now useful to identify the form of the perturbations x and ^+ that are 
pure gauge transformations of the string configuration. It is easy to see that per- 
turbations of the form 

5cP - Z(7V0O, dW, = -iDo,^P , (3.29) 

where ip is a real L(SU(2)) valued function and the subscript denotes the un- 
perturbed fields, is an infinitesimal gauge transformation of the original vortex 
solution. (In (3.29), Wi represents r • Wj.) If we now require that these purely 
gauge perturbations do not afi^ect the string configuration itself, then we can only 
have 

'/'=''('-)(_i% %**) (3-30) 
where s{r) is any smooth function. This means that perturbations given by 

x \ _ ^(^\ ( -gf 



are pure gauge perturbations that do not affect the string configuration. Therefore, 
such perturbations cannot contribute to the energy variation and must be annihi- 
lated by O. Then, in the two-dimensional space of (X)^+) perturbations, we can 
choose a basis in which one direction is pure gauge and is given by (3.31) and the 
other orthogonal direction is the direction of physical perturbations. The physical 
mode is, 

C = (1 - iv)x + • (3.32) 

It was a good check on our algebra that on eliminating ^-|- in terms of ( and x 
in eq. (3.28) the functional reduces to one depending only on (: 

SE[C] = 2tt j dr r(OC (3.33) 

where O is the differential operator 

«-J|:(i^)--« 
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and 

Uir) = + + j • (3-35) 

To summarise, the question of stability reduces to asking if the operator O 
has negative eigenvalues in its spectrum. That is, whether the eigenvalue u of the 
Schrodinger equation, 

OC = a;C, (3.36) 

can be negative. The eigenfunction ( must also satisfy the boundary conditions 
C(r = 0) = 1 and C — c (c is some constant) as r — > oo. 

In this way we have reduced the stability analysis down to one Schrodinger 
equation which we will solve numerically. 

4. Numerical Analysis 

We solve eq. (3.36) together with the Nielson-Olesen eqs. (2.9) and (2.10). How- 
ever, it is convenient to work with rescaled dimensionless variables. Hence, we 
define 

P^^f, V^%, R^^r. (4.1) 
V 2 ' 2^2 ^ ^ 

In terms of these dimensionless variables, the Nielsen- Olesen equations (2.9) and 
(2.10) become, 

V" + 2(1 - V)P^ = (4.3) 
R 

where primes now denote differentiation with respect to R. The functions P and 
V also satisfy the boundary conditions: 

P(0) = = F(0), P(oo) = 1, V{oo) = 1 (4.4) 
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The problem now has only two free parameters: cos 9^ and This may be seen by 
rescaling fields and coordinates in the operator O as in (4.1). With these rescalings 
the quantities 8+ and P+ in the eigenvalue problem (3.36) get replaced by: 



The rescaled eigenvalue problem was solved by using a fifth order Runge-Kutta 
algorithm. We kept /3 fixed and found 9w for which the lowest eigenvalue changes 
sign. We repeated this procedure for several values of /3 and found the correspond- 
ing values of critical parameters (y^, sin^ 6*^;). The above method was used to 
scan the range 0.07 < /? < 1.0. Lower values of P make the numerical analysis 
fairly intensive since then there are two widely different scales in the problem cor- 
responding to the two widely different masses. Our results are shown in Fig. 1 
where we plot the critical values of (the ratio of the Higgs mass to the Z mass) 
versus the corresponding values of sin^ 9w In sector III, on the right-hand side of 
the data line, equation (3.36) had no negative eigenvalues implying string stability. 
Thus we may distinguish three sectors in Fig. 1: sector I where the electroweak 
strings are unstable, sector III where strings are stable, and, the presently unex- 
plored region shown as sector II (/? < 0.07 or mjj < 2AGeV). It is evident that 
the physically realized values: sit? 9^ = 0.23 and = > 0-62 (see Ref. 

14) lie entirely inside sector I. This brings us to the main result of this paper: if 
the standard electroweak model is the physically realized model, then the existing 
vortex solutions in the bare model are unstable. 




where. 



P; = (1 - 2cos'^9wVf + 2cos^9wB?P'^ 
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5. Outlook 

We have rigorously established that the electroweak model admits stable vortex 
solutions for a certain range of parameters. This result is exciting in that it makes 
the possibility of observing coherent states in particle physics closer to reality. On 
the other hand, it is somewhat dissappointing that the values of the parameters that 
Nature has actually chosen are such that the vortex is unstable. However, this still 
does not mean that the vortex will be unstable in the real world since our analysis 
only applies to the bare electroweak model. In the context of the early universe, 
for example, we must do a stability analysis at high temperatures. One should 
also consider the possibility that Nature has chosen an extension of the standard 
model where the Higgs potential is more complicated. In such circumstances, the 
stability issue would have to be readdressed. 
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FIGURE CAPTIONS 



A map of parameter space showing the results of the stabihty analysis. Sector I 
contains unstable strings, sector III contains stable strings and we have not explored 
sector II. We also indicate the physically allowed range of parameter space. The 
data from LEP constrains the Higgs mass to Mjj > bSGeV which implies > 
0.62GeV and the observed sin^ is 0.23. 
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